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Abstract- This study presents the technology to combine the 
remote sensing image of SPOT and FORMOSAT-2 satellite 
image by Fuzzy theory and Back Propagation Network (BPN). 
This method adopt three experience identify factors of NDVI, 
shape, and color to establish the membership function. The fuzzy 
results show that the successful rate of identification was about 
87 percent. The BPN results show that the successful rate of 
identification was about 91 percent. Therefore the result proves 
that the three indexes are the best choice in extracting of change 
area from image, because the error is the lowest among the two 
methods in this study. 

Keywords- Component; Fuzzy Theory; Image Extraction; 
Image Variation; Image Identification; Back Propagation Network 
(BPN) 

I. INTRODUCTION 

The objective of monitoring the change of catchment area 
is to avoid damage from nature disaster or excessive 
development. Previous method uses software to identify the 
build land category, followed by human experience to check it 
one by one. However, it needs a lot of time to identify image 
from remote sensing images. If it is possible to automatically 
study and train the capability of selection system, then the 
time of distinguish can be effectively reduced. 

For example, Meng-Lung Lin and Cheng- Wu Chen used 
fuzzy model and remote sensing data about Normalized 
Difference Vegetation Index (NDVI) in the winter, spring, 
summer, and fall of 2000 and 2005. Their research analyzes 
land cover maps and landscape sensitivity, and then estimate 
the sensitivity of the ecosystems in the landscape [1]. Natalie 
Campos etc., 2011 used High-resolution imagery and find it 
were successfully assess map changes in land-cover patterns in 
research area. However, previous these researches in similar 
area were comparing lacking. They used NDVI image 
differencing and Classification Tree (CT) methods in their 
study. Then analyze to connect time with the changes in 
spatial distribution of mountain resort development (MRD)[2]. 

Using image segmentation software to produce image 
segmentation vector and extract the variation region is a 
common practice for remote sensing technique in the past. The 
process is: selecting image to segment, outputting the 
segmentation vector file, using the software to extract the 
vector layer, and screening for variation area. The work is 
done by the accumulation of human experience. The memory 
capability of human can accumulate the previous experience, 
and it can compile the things touched to the rules of 
information. 

Human brain can trigger the corresponding rules to the 
next move. If there were no corresponding rules, it can also 
trigger the fuzzy space for self-learning and reasoning. This 



capability can be simulated by the case base reasoning with 
fuzzy theory. Therefore, we use the fuzzy theory and case base 
reasoning to construct a quantitative fuzzy space structure to 
learn the human brain logic and estimate the variation region 
in the remote sensing image. 

II. RESEARCH AREA 

The research area locates in Shihmen Reservoir in 
Taoyuan County (Fig. 1). Base on the report (Level 3 
Management of Shihmen Reservoir) from Water Resources 
Agency, Ministry of Economic Affairs, it can be seen that the 
areas of orchard expand rapidly and mainly spread on sloping 
fields at altitude of 800 meters which is good for planting 
temperate fruit trees (peaches) but lower altitude areas are for 
tangerines. The usage of lands follows by the demand of 
tourism and economy benefits. It has been a significant change 
compared with earlier types of crops; therefore it is essential to 
monitor the environment of Shihmen Reservoir. 




Fig. 1 Shihmen Reservoir 



If telemetry can be used to monitor this area, it can quickly 
browse the full view of reservoir lands and strengthen the 
ability of monitoring variation region instead of spending great 
amount of time and money by using manpower. It can 
improve the traditional way of monitoring to gain efficiency. 
What's more, there have drafted down some preventive 
measures of man-made or nature disaster to reduce the 
influence from disaster and ensure lands can last forever. 

III. METHODOLOGY 

Three indexes were adopted in this study as experience 
rules for quick discrimination, they were: differences in NDVI 
values, shape ratio, and color difference values (Table I and 
Figs. 2-4). 
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TABLE I MEMBERSHIP FUNCTION AND NUMERICAL VALUE OF THE THREE 
EVALUATION INDEXES 



Evaluation 
Index of 

Variation 
Regional 


Type if 
Membership 
Function 


Attribute if 
Fuzzy 
Interval 


Membership 

Function 
Value 


Differences 
in NDVI 
Values 




NDVI 

Differences 
Values is 
Positive 


0.3-2 


Trapezoidal 


NDVI No 
Change 


-0.3 — 0.3 






NDVI 
Differences 
Values is 
Negative 


-2--0.3 






No Change 


1 To 1.5 
Times 


Shape Ratio 


Trapezoidal 


Diffusion 
Change 


1.5 Times To 
10 Times 






Reduction 
Change 


0—1 Time 






Similar 


-30-30 


Color 
Difference 
Value 


Trapezoidal 


Dissimilar, 
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Dissimilar, 
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Fig. 2 Differences in NDVI values 



Ret uuudmnge 



0.8 

a 

Z 0.4 

0.2 



Diffusion 



4 6 
shape 
Fig. 3 Shape ratio 



Hi 



0.8 



S3 

■a o.6 



0.1 



0.2 



Dissimilar 



Similar 



Dissimilar 



— i 1 — 



— I L_ 



-600 -400 -200 200 400 600 
color 

Fig. 4 Color difference values 

A. Differences in NDVI Values 

Chao-Yuan Lin etc., 2010 used multi -temporal satellite 
images index (NDVI) to automatically extract landslide and 
image classification. Usually, image classification uses NDVI 
to increase accuracy. Because it can effective to classify the 
bare land and vegetation cover [3]. The a value, NDVI [4], is 
to use the characteristic of green plant which can absorb the 
red light and reflect the near red wave. According to wave 
length from the remote sensing data, it can provide the plant 
growth trend for the variation area. The meaning of NDVI 
value difference is to compare the ratio of red light absorption 
and near-infrared light reflection by green plant in the remote 
sensing image. The numerical value of NDVI difference is the 
difference of the early image and late image. Its value was 
ranged between -2 — 2. 

1 ) Define the NDVI value of previous image is a, and 
the value of post-image is b. The range of a single NDVI value 
is between -1-1. 

2) Let f(a) = {1,-1} f(b)= {-1,1} .then f(b)-f(a) = -2 
and 2 , so the maximum threshold value is between -2-2. 

3) The different value of NDVI is in the range between 
-2-2. According to the induction of manual experience, if the 
difference value of NDVI is lower than 0.3, it could be 
regarded as having no change. 

4) Let f(a)= {0.2,0.5} f(b)= {0.5,0.2} , we set the 
range as if the value is between -0.3 - 0.3, then there is no 
change. If the value is between 0.2-0.3, then there is change. 

The function of NDVI is to rapidly analyze the plant 
growth rate in the variation region. Otsu has proposed a rule 
for plant in the remote sensing image in 1979[5].The color 
value of image was divided into two areas as vegetation and 
non-vegetation. 



NDVI 



(NIR - R) 

(nir + R) 



(1) 

Where, NIR is the near infrared band and R is the red band. 

B. Shape Ratio 

The shape ratio (a P) is the area after the event divided by 
the area before the event. The outer shape is used to analyze 
whether it has achieved the variation, as shown in Fig. 5. 
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Old Image (2008.07.21) New Image (2008.09.21) 




Fig. 5 Shape ratio 

The definition of shape ratio is the ratio of post stage area 
to the previous stage area. It has to be assigned with the value 
according to the error factors of location, shape, and twist, etc. 
If the value is lower than 1, it is the reduction change. If the 
value is 1 to 1.5, then there is no change. If the value is greater 
than 1.5, then it is the diffusion change. 

The shape ratio is not applicable to the land use type as: 
from never to have, from having to no, and the simplest 
change of type. The shape ratio is applicable to the condition 
when the land use types have no change, and only the area 
shape changes. The shape ratio is also applicable to the area 
which has been labeled as the variable region. 



Shape ■■ 



ifi-a) 
to 



(2) 



C. 



Color Difference Values 

The color difference is the color before and after the event. 
In any pixel of the remote sensing data, the maximum value of 
RGB will not exceed 255+255+255 = 765, as shown in Fig. 6. 
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Fig. 6 Color difference values 
IV. CONSTRUCT FUZZY 

The MATLAB fuzzy toolbox was used to construct a 
structure of rapid extraction. Three groups of indicators were 
identified. The functional forms and membership function 
value were shown in Figs. 4-7. Since every index has three 
attribute of fuzzy interval. So every index has four critical 
ranges. Kaufmann and Gupta [6] induced several types of 
membership function. Most suitable for all indexes is the 
trapezoidal fuzzy number defined as: 



f(x,a,b,c,d) = 






x < a 


x — a 






a < x <b 


b — a 




I 1 


b< x <c 


d — x 


c<x<d 




d — c 







x > d 



a < b,c < d 

The trapezoidal function was chosen in this study because 
the interval of threshold value for the three indexes in this 
study is relatively fixed. Twenty seven rules of < IF... THEN > 
were formulated as the following: 

1) If (NDVI is no change) and (shape is no change) and 
(color is Similar) then (expected value is no change) 

2) If (NDVI is no change) and (shape is no change) and 
(color is Dissimilar) then (expected value is yes change) 

3) If (NDVI is no change) and (shape is no change) and 
(color is Dissimilar) then (expected value is yes change) 

4) If (NDVI is no change) and (shape is Reduction) and 
(color is Similar) then (expected value is yes change) 

5) If (NDVI is no change) and (shape is Reduction) and 
(color is Dissimilar) then (expected value is yes change) 

6) If (NDVI is no change) and (shape is Reduction) and 
(color is Dissimilar) then (expected value is yes change) 

7) If (NDVI is no change) and (shape is Diffusion) and 
(color is Similar) then (expected value is yes change) 

8) If (NDVI is no change) and (shape is Diffusion) and 
(color is Dissimilar) then (expected value is yes change) 

9) If (NDVI is no change) and (shape is Diffusion) and 
(color is Dissimilar) then (expected value is yes change) 

10) If (NDVI is negative) and (shape is no change) and 
(color is Similar) then (expected value is yes change) 

11) If (NDVI is negative) and (shape is no change) and 
(color is Dissimilar) then (expected value is yes change) 

12) If (NDVI is negative) and (shape is no change) and 
(color is Dissimilar) then (expected value is yes change) 

13) If (NDVI is negative) and (shape is Reduction) and 
(color is Similar) then (expected value is yes change) 

14) If (NDVI is negative) and (shape is Reduction) and 
(color is Dissimilar) then (expected value is yes change) 

15) If (NDVI is negative) and (shape is Reduction) and 
(color is Dissimilar) then (expected value is yes change) 

16) If (NDVI is negative) and (shape is Diffusion) and 
(color is Similar) then (expected value is yes change) 

17) If (NDVI is negative) and (shape is Diffusion) and 
(color is Dissimilar) then (expected value is yes change) 

18) If (NDVI is negative) and (shape is Diffusion) and 
(color is Dissimilar) then (expected value is yes change) 

19) If (NDVI is positive) and (shape is no change) and 
(color is Similar) then (expected value is yes change) 

20) If (NDVI is positive) and (shape is no change) and 
(color is Dissimilar) then (expected value is yes change) 



URSA Vol.2 Iss. 2 2012 PP.30-35 www.ijrsa.org ©World Academic Publishing 

-32- 



International Journal of Remote Sensing Applications 



URSA 



21) If (NDVI is positive) and (shape is no change) and 
(color is Dissimilar) then (expected value is yes change) 

22) If (NDVI is positive) and (shape is Reduction) and 
(color is Similar) then (expected value is yes change) 

23) If (NDVI is positive) and (shape is Reduction) and 
(color is Dissimilar) then (expected value is yes change) 

24) If (NDVI is positive) and (shape is Reduction) and 
(color is Dissimilar) then (expected value is yes change) 

25) If (NDVI is positive) and (shape is Diffusion) and 
(color is Similar) then (expected value is yes change) 

26) If (NDVI is positive) and (shape is Diffusion) and 
(color is Dissimilar) then (expected value is yes change) 

27) If (NDVI is positive) and (shape is Diffusion) and 
(color is Dissimilar) then (expected value is yes change) 

The fourth stage investigation data from the project of 
SWCB-97-157 in 2008 were taken as the experience cases. 
Ho-Wen Chen etc., 2009 developed a neural-fuzzy inference 
approach with used the FORMOSAT-2 (8m spatial resolution) 
multi-spectral satellite image. It can to analysis the land use 
and land cover patterns in a fast growing city. Their each 
fuzzy-neural rule incorporated with both the satellite image 
spectral and textural features. So, it increased the advantages 
of classification efficiency and accuracy [7]. 

However, in this paper try to use better high resolution 
satellite fusion images (FORMOSAT-2 2m spatial resolution 
and SPOT) with three indexes, the results were further verified. 
The expected value of variations was calculated by the module 
and obtains the value of 0.67. The module chart is shown in 
Fig. 7. It means that if there is a possibility of change, its value 
will be greater than 0.67. Otherwise, the value of change will 
be lower than 0.67. The results of experience cases were 
shown in Table II, III. From the table results, there are three 
cases with the expected value of variance lower than 0.67 
among the 32 cases. So the success rate of distinguish is (55- 
7)/55=0.87. The recognition rate is about 87%. 



TABLE II THE THIRD EXPERIENCE CASE IN 2008 
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Fig. 7 Module chart 
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TABLE III THE FOURTH EXPERIENCE CASE IN 2008 



No 


Differences in 
NDVI Values 


Shape 
Ratio 


Color 
Difference 
Value 


Expected 
Value of 
Variations 


Dl 


-0.25 


1.00 


-128.00 


0.67 


D2 


-0.63 


1.00 
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0.67 


D3 
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0.67 


D4 
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V. NEURAL NETWORK MODEL 

One neural network model, the Back Propagation Network 
(BPN) was used in this study to find the best results [8] . The 
experimentation data of BPN are to the real value (Table IV). 
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TABLE IV THE THIRD AND FOURTH ERROR RANGE 





Uitterences 

in INDV1 

Values 


Shape 
Ratio 


Color 
Difference 
Values 


Assessment 
Index 


Expected 
Value 


Error 
Range 

\ fO) 


CI 


-0 29 


: 


89 


19 73 


19 51 


1 14% 


CO 


-u.zo 





i ftn 

lOU 


OA 10 
Z4-. 1Z 


93 so 

Z J. JZ 


o sn^i, 

Z.JU /o 


P3 


n i s 

-U. 1 J 


: 


ftO 

-QZ 


1 A Al 
14-. 4- / 


1 S 43 
1 J.4- J 


ft ftft^ 

o.oo /o 


CA 







83 

-0 J 


00 84 
ZZ. O^f 


i y. jo 


1 S 04"% 


L.J 


-U.4Z 


: 


1 01 


OS 47 
ZJ.4 / 


94 48 
Z4.4fi 


1 Oftc,' 
J.fiO /o 


rv, 
L.O 


n n i 


- 


OJ 


11.1/ 


1 1 QA 
1 1.84 


^ mo/ 
o.Ujvo 


CI 


n no 

-U.Uo 





o 


y.uo 


1 O £0 
IZ.OZ 


io ino/ 
jy.jUvb 


r^o 
L.5 


n no 
-u.uy 





1 11 
1 JJ 


1 /.Ol 


17 10 
1 /. 1Z 


O 11C/ 
L.I 170 


fA 


-U.UO 


J.OO 


-ZJ 


4U.8 / 


1/1 07 

j4.y / 


1/1 /I /I 0/ 

14.44vo 




A A'C 
U.Uj 


1.60 


-o 


1 A A/1 

14.U4 


1 A £.C\ 

14.0U 


1 A CO/ 
J.y JVC 


pi i 

L.1 1 


n no 
U.Uo 




00 
-ZZ 


1U. 1U 


1 O Q1 
1Z.8J 


O/^ QOO/ 

zo.yyvo 


L.1Z 


n i o 
u.iy 


— - — 


-yo 


1 1 

i /.oy 


1 Q 

18.80 


O.OZvo 


L.1 J 


n i 
U.j 


— - — 


1 1 1 
-in 


1 /10 
Z1.4Z 


on 

ZU.fiJ 


O £SO/ 

Z.ujvo 


L.14 


n m 
U.U / 


— - — 


OS 
-OJ 


i j.y / 


1 7 oy 

1 /.Zo 


01 A7C/ 

Zj.o/ To 


pic 

L.1 j 


a n 
U.z / 


— - — 
— J — 


OA 

zy 


1 J.jl 


i a 


1 A "70 0/ 

jU. /fi% 


pi/; 


n os 

U.Zj 




n i 
i j 


1 1 £Q 

i /.oy 


O 70 

y. /y 


/I /I £ S 0/ 
44.0 JVC 


C17 


A 1 £ 

U.lo 


— - — 
— J — 


A C 

4j 


1 J.Ol 


1 A /I A 

lU.4o 


Ol 1 1 07 
Z J .11% 


PIQ 

L.18 


n i o 

U.lo 




1 1 
1 J 


i o no 
1Z.UZ 


1 1 01 
1 1.Z1 


0. / Zvo 


PIQ 

L.iy 


n no 
U.Uo 


Q7 

z.y / 


00 


0£ 1 1 

Zo. 1 1 


OS QQ 
Zj.fio 


n ono/ 
U.yUvo 


pin 


0.25 


— J — 


-53 


16.38 


15.50 


5.39% 


P1 1 
L.zl 


n i s 

U.l 3 




hi 

-ID 


1 < 1 Q 

i j. iy 


16.65 


Q ££0/ 

y.oovo 


poo 

L.ZZ 


-0. 17 


— - — 
— J — 


171 


22.09 


21.80 


1 .33% 


r\i 
JJ1 


n os 

-U.Z J 




1 00 
-IZo 


O 1 OC 
Zl.Zfi 


1 O 07 

iy.fi / 


0.0470 


TOO 

DA 


n ^i 
-U.Oj 


— - — 


-5 
-J 


OO ^1 
ZZ. Ol 


01 £Q 

Zj.oy 


/I lie/ 
4. / / To 


n3 


-0 62 


— - — 


-13 


23 02 


23 22 


U.7U /O 


D4 


-0.47 


: 


-22 


19.85 


19.31 


2 73% 


ns 


-0 22 


: 


Q3 

7J 


18 25 


17 53 


3 Q0% 


D6 


-0.34 


: 


-53 


18.63 


16 75 


1 08% 


D7 


-0.5 


: 


_\ 


19.23 


20.59 


1 07% 


D8 


-0 38 


2 56 


-26 


28 27 


08 S4 

ZO. JH- 


n Qft% 

U.7U /O 


D9 


-0.14 




-90 


16.05 


17 19 


7 12% 


U1U 


f\ AO 
-U.4Z 


— - — 


is 
-jj 


i y.4j 


i s no 

1 fi.UZ 


7 IStK, 
/.JJ /o 


Dl 1 


n 1 
-U.l 


— - — 


Al 
-4 / 


1 O O/l 
1Z.Z4 


1 A A ^ 
14. 4J 


1 q n^o/ 
Ifi.Uovb 


D1Z 


n /i 
-U.4 


— - — 


OQ 

zy 


1 Q 

Ifi.JO 


iy. lo 


1 O/l 07 
J.Z470 


D13 


-U.JJ 


— - — 
— J — 


OO 


19. /j 


i n to 
19. /z 


U.U4vo 


TO 1 A 
D14 


-U.jl 




1 C\ 

1U 


on m 
ZU.U / 


1 01 

Z1.Z1 


s 7no/ 
j . / U/o 


Dl J 


A C/l 

-U.j4 


— - — 
— J — 


-16 


O 1 O 1 

zl.zl 


O 1 1 A 

zi.iy 


n 110/ 
U.l jye 


DIO 


n /10 
-U.4Z 




1 A 
-14 


i o nc 
Ifi.Ufi 


1 07 
Ifi.Z / 


1 noo/ 
1 .UZvo 


nn 
Dl 1 


-U.Jo 


— - — 


10 
-JO 


01 

Zj.oj 


O 1 77 


7 O/l O/ 
/ .y47o 


D18 


A /I A 

-U.44 


— - — 
— J — 


Ol 
-ZJ 


19. 17 


1 O CO 

Ifi.Jfi 


1 mo/ 
j.U/vo 


TO 1 O 
Diy 


n i 




1 
-1Z 


1 A 

14. yj 


1 S 7C 
1 J. /J 


s 100/ 

J.JzTc 


TO O A 

DzU 


A A A 

-U.44 


— - — 
— J — 


OA 

2U 


18.97 


in oi 
19. ol 


/I /I IO/ 

4.4 j 70 


r*o 1 
DZ1 


Ot A 1 

-U.41 




i nn 
-1UU 


01 A*Z 

Zj.4j 


1 o oc 

ly.zj 


1 7 ono/ 

1 / .yuvo 


DZZ 


n ac 
-U. jj 


— - — 


1 j 


1 ^ /in 
10.4U 


1 7 17 
1 /. J / 


s 1 0/ 
j.y 1 To 


D23 


-0 27 


: 


Q0 

- "Z 


19 43 


17 68 


9 00% 


D24 


-0 22 


: 


-62 


16 22 


15 80 


2 59% 


D25 


-0 42 





41 


1 V.O J 


on 03 

ZU.Z J 


1 Q3% 

1 .V J /o 


D26 


n sft 

-W.JO 


: 


-75 


OS S7 
ZJ. J / 


21 21 


1 7 ns% 

1 / .OJ 10 


D27 


-0.56 




-40 


23.28 


21.25 


8.71% 


D28 


-0.43 




-20 


18.72 


18.39 


1.76% 


D29 


-0.45 




-47 


20.99 


18.65 


11.14% 


D30 


-0.28 




-19 


14.91 


15.39 


3.20% 


D31 


-0.3 




26 


15.87 


16.55 


4.34% 


D32 


-0.57 




-3 


21.11 


22.28 


5.53% 


D33 


-0.05 




-104 


14.71 


18.36 


24.80% 


Mean 


8.96% 



A. Training and Test Model 

In this mode, the original data was used to establish the 
model. The input was the unit for training, the number of 
hidden layer, learning cycle, and speed of learning cycle, etc. 

B. Authentication Model 

The authentication mode is to use the non-training data to 
evaluate the reliability of the model. The results can be the 
information of the scattering plot, matrix, or error analysis. 



The output variable will also be estimated to compare their 
error. 

C. Inference Model 

The inference model is used to predict certain condition 
that we want to know. The experienced value was derived 
from the trained mode. (citation : 

http://www.wisdomsoft.com.tw) 

The value of assessment index is calculated as the 
following: 

. , [Differences in NDVI values) ,„„ 
assessmentindex = '- x 100 + 



(Shape ratio) (Color difference value) 

-xl00 + - -xlOO 



15 



1530 



(3) 



The prediction results of these BPN neural network 
schemes were shown in (Table IV). The error of BPN to the 
real value is 91.04%. 

VI. CONCLUSION 

It is necessary to appropriately monitor the environment 
according to the characteristic of land use in the watershed. 
Remote sensing technology helps to provide the real time data, 
and can quickly glimpse the land picture status of watershed. 

This technology enforces the capability of variation region 
monitoring, and could save time and cost. The result reveals 
that the FUZZY and BPN have better feasibility of predictive 
values, which proves that the three index is better choice in 
extracting change area from image. Therefore, we will use this 
experience to identify factors of NDVI, shape, and color to 
establish the membership function basic module for further 
research. 
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